Scattering theory of topological invariants in nodal superconductors 
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Time-reversal invariant superconductors having nodes of vanishing excitation gap support zero- 
energy boundary states with topological protection. Existing expressions for the topological invari- 
ant are given in terms of the Hamiltonian of an infinite system. We give an alternative formulation 
in terms of the Andreev reflection matrix of a normal-metal-superconductor interface. This allows 
to relate the topological invariant to the angle-resolved Andreev conductance, also when the bound- 
ary state in the superconductor has merged with the continuum of states in the normal metal. A 
variety of symmetry classes is obtained, depending on additional unitary symmetries of the reflec- 
tion matrix. We derive conditions for the quantization of the conductance in each symmetry class 
and test these on a model for a 2D or 3D superconductor with spin-singlet and spin-triplet pairing, 
mixed by Rashba spin-orbit interaction. 
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I. INTRODUCTION 

The topological classification of superconductors relies 
on the existence of an excitation gap in the bulk of the 
material, that prevents transitions between topologically 
distinct phases. 1 ' 2 The gap of a topological superconduc- 
tor closes only at the boundary, where propagating states 
with a linear dispersion appear. The protected boundary 
states are counted by a topological invariant Q, expressed 
either in terms of the Hamiltonian of an infinite system 3 
or in terms of the scattering matrix for Andreev reflection 
from the boundary with a normal metal. 

Nodal superconductors with time-reversal symmetry 
also have boundary states, forming flat bands in the 
middle of the bulk gap. 5 The same topological consid- 
erations do not apply because the gap vanishes in the 
bulk for certain momenta k on the Fermi surface (nodal 
points). Examples include the cuprate superconductors 
(gap oc k x k y )® and a variety of superconductors without 
inversion symmetry. 7 Nodal superconductors may also 
appear as an intermediate phase in the transition from a 
topological superconductor to a trivial one. 8 ' 9 

A topological invariant can still be constructed in 
a nodal superconductor for a translationally invariant 
boundary, 10,11 conserving the parallel momentum fen. 
The value of Q(feii) can only change if fen crosses a nodal 
point. This topological invariant again counts the bound- 
ary states, which are now non-propagating dispersionless 
states (pinned to E = for a range of feii). 

In Refs. 10,11 the topological invariant Q(feii) of a 
nodal superconductor takes the form of a winding num- 
ber, calculated from the Hamiltonian of a translationally 
invariant infinite system. Here we present an alternative 
scattering formulation, which expresses Q(feii) as a trace 
of the Andreev reflection matrix. Since the conductance 
of a normal-metal-superconductor (NS) interface is ex- 
pressed in terms of the same Andreev reflection matrix, 
this alternative formulation allows for a direct connec- 
tion between the topological invariant and a transport 



property. 

If the NS interface contains a tunnel barrier, the 
angle-resolved conductance G(ku) measures the density 
of states and directly probes the flat surface bands as a 
zero-bias peak. 12 For a transparent interface the bound- 
ary states in the superconductor merge with the contin- 
uum in the metal, resulting in a featureless density of 
states, but the zero-bias peak remains. 13 Here we relate 
the height of this zero-bias peak to the value of the topo- 
logical invariant. While in general this relation takes the 
form of an inequality, a quantized conductance, 



G(fe,|) = |Q(fe|,)| x 2e 2 /h, 



(1.1) 



may result under certain conditions which we identify. 

The outline of this paper is as follows. In the next sec- 
tion we formulate the scattering problem and construct 
the topological invariant from the Andreev reflection ma- 
trix. We make contact in Sec. Ill with the Hamiltonian 
formulation, by closing the system and showing that we 
recover the number of flat bands at the boundary. We 
then return to the open system and in Sec. IV relate the 
angle-resolved zero-bias conductance to the topological 
invariant. So far we only assumed the basic symmetries of 
time-reversal and charge-conjugation. The effects of ad- 
ditional unitary symmetries are considered in Sec. V. We 
apply the general theory to a model of a two-dimensional 
(2D) nodal superconductor in Sees. VI and VII, includ- 
ing also the effects of disorder. Effects that are specific 
to 3D are discussed in Sec. VIII. We conclude in Sec. IX. 



II. TOPOLOGICAL INVARIANT FOR 
ANDREEV REFLECTION 

A. Chiral symmetry 

We study the Andreev reflection of electrons and holes 
at the Fermi level from a planar interface between a nor- 
mal metal (N) and a superconductor (S). (See Fig. 1.) 
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FIG. 1: Interface between a superconductor (S) and a normal 
metal (N). The reflection matrix r(feii) relates the amplitudes 
of the incident and reflected waves (arrows, both normal re- 
flection and Andreev reflection are indicated). The conduc- 
tance of the NS interface is measured by applying a voltage 
difference V between the normal metal and the grounded su- 
perconductor. 



The component fen along the interface of the momentum 
fc is conserved, so we can consider each fen separately 
and work with a one-dimensional (ID) reflection matrix 
r(fc||). For fc not in a nodal direction (nonzero excitation 
gap) this is a unitary matrix, 



r(feii)r' (feii) = 1. 



(2.1) 



The dimension of the reflection matrix is 4 x 4, with 
basis states (ip e t,ipei,iph^,iphi) labeled by the spin f,4- 
and the electron-hole e, h degrees of freedom. The e, h 
grading produces four 2x2 submatrices, 

Normal reflection (from electron to electron or from hole 
to hole) is described by r ee and rhh, while The and r e h 
describe Andreev reflection (from electron to hole or the 
other way around). 

The two fundamental symmetries that we impose are 
time-reversal and charge-conjugation symmetry. Time- 
reversal symmetry requires 



r(fe||) = cr a r T (-fc||)cr y , 



(2.3) 



while charge-conjugation symmetry at the Fermi level re- 
quires 



r(fc||) = ^r* (— fel^Taj 



(2.4) 



The Pauli matrices Ui and r, act on, respectively, the spin 
and electron-hole degrees of freedom. (For later use we 
denote the 2x2 unit matrices by <7o and tq.) 

Taken together, Eqs. (2.3) and (2.4) represent the chi- 
ral symmetry relation 



r(fe||) = (a y <E> r x )r ] (k\\)(a y <X> t x ). 



(2.5) 



This is the ID symmetry class AIII in the periodic table 
of topological phases. 3 



It is convenient to represent the symmetry relations in 
terms of the matrix -R(feii) = {a v ® r x )r(k\\), which is 
both Hermitian and unitary, 



R = R\ R 2 = 1. 



The submatrices in Eq. (2.2) appear in R as 
i?(fc||) 



Rhe(k\l) Rhh{k\\ 
Ree(k\\) Reh{k\\ 



(2.6) 



(2.7) 



where R pq = <J y r pq . The two blocks i?/ te and R e h are 
Hermitian, while R ee = R 



■hh- 



B. Topological invariant 

The Z topological invariant of ID reflection matrices 
in class AIII is given by 14,15 



Q(fe,|) = iTrfl(fc,|) 



■ Tr cr y [r fce (fc||) + r eft (fc||)]. 



(2.8) 



In view of Eq. (2.6), the 4x4 matrix R has eigenvalues 
±1, so the value of Q G {-2,-1,0,1,2}. This value is 
fe||-independent as long as the reflection matrix remains 
unitary. For fc in a nodal direction, the reflection matrix 
is sub-unitary and the topological invariant may change. 
Application of Eq. (2.3) gives the relation 



i?(-fc,|) 
which implies that 



-T x R T (ku)T x , 



G(-fcii) = -C(fcii) 



(2.9) 



(2.10) 



If feii =0 one necessarily has Q = 0. For this time- 
reversally invariant momentum the Pfaffian of the anti- 
symmetric matrix a y r(0) (equal to ±1) produces a Z2 
topological invariant, 14,15 characteristic of the ID sym- 
metry class DHL We write this invariant in the form 



Qo = 1 + Pf<v(0) g {0,2}, 



(2.11) 



so that for Qo, as well as for Q, the value indicates the 
topologically trivial phase. 



III. TOPOLOGICALLY PROTECTED 
BOUNDARY STATES 

The scattering formulation of topological invariants 
refers to an open system, without bound states. In the 
alternative Hamiltonian formulation, the topological in- 
variant counts the number of dispersionless boundary 
states (flat bands at the Fermi level, consisting of edge 
states in 2D or surface states in 3D). 10 ' 11 ' 16-19 To relate 
the two formulations, we close the system by means of 
an insulating barrier at the NS interface, and show that 
|<2(fc||)| boundary states appear. 
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The calculation closely follows Ref. 14. The number of 
boundary states at fey equals the number of independent 
solutions ip of 

[l-n(fc||)r(fc||)]V = 0. (3.1) 

The unitary matrix r\ is the reflection matrix of the bar- 
rier, approached from the side of the superconductor. We 
can write this equation in terms of Hcrmitian and unitary 
matrices R\ — r\(a y ® t x ) and R 2 = {cr y <8> T x )r, which 
we decompose as 

Rl = U iDi ul A = ( la J* ^ J . (3.2) 

(The notation 1 m indicates the M x M unit matrix and 
U\, U 2 are unitary matrices.) Eq. (3.1) takes the form 

(1 - DiUD 2 U^)i/j' = 0, (3.3) 

with U = U\U 2 and V' = E/fy. 

We decompose U into N x M submatrices A N>M , 

jj _ (M+Q 1 ,2+Q 2 A 2 +Q 1 ,2-Q 2 \ /g 4 n 

\^2-Qi,2+Q 2 ^2-Q 1 ,2-Q 2 /' 

Since 

U-D 1 UD 2 = 2(. ° ^+q,,2-q 2 A (35) 
we can rewrite Eq. (3.3) as 

\a 2 -q u2 +q 2 o y 
with ip" = u\i). 

For any matrix Ajv,m with N < M there exist at 
least M — N independent vectors v of rank M such that 
A NM v = 0. Therefore Eq. (3.6) has at least \Qi + Q 2 \ 
independent solutions. These are the topologically pro- 
tected boundary states. 

Because the insulating barrier is topologically trivial, 
Qi = 0, while Q 2 = Q is the topological invariant of 
the superconductor, so it all works out as expected: The 
topological invariant of the open system counts the num- 
ber of boundary states that would appear if we would 
close it. 

Both values Q and — Q of the topological invariant 
produce the same number Af = \Q\ of boundary states 
if the superconductor is terminated by a topologically 
trivial barrier (an insulator or vacuum). The sign of 
the topological invariant matters if wc consider the in- 
terface between two topologically nontrivial supercon- 
ductors 1,2. The combined number of boundary states 
-A/total = I Qi + Q.1 1 = |-A/"i ± M 2 1 is the sum or difference of 
the individual numbers depending on whether the topo- 
logical invariants have the same or opposite sign. 



IV. RELATION BETWEEN CONDUCTANCE 
AND TOPOLOGICAL INVARIANT 

By considering an open system when formulating the 
topological invariant, we can make direct contact to 
transport properties. The angle-resolved zero-bias con- 
ductance of the NS interface is given by 

G(fe||) = G Trr, e (fe||)rL(fe||), (4.1) 

with Go = 2e 2 /h the Andreev conductance quantum. We 
wish to relate this transport property to the topological 
invariant (2.8). 

For that purpose it is convenient to work with the ma- 
trices Rhe = cryrhe and R e h = <y y r e h 1 since these are Her- 
mitian (unlike the r^e and r e h themselves). For brevity 
we omit the label fe||. The squares R\ e and R? eh have the 
same set of Andreev reflection eigenvalues p n € [0,1], 
which are also the eigenvalues of r he rj le . 

On the one hand we have the conductance 

G/G = TrR 2 he = TrR 2 eh , (4.2) 

and on the other hand the topological invariant 

Q = \TT(R he + R eh ). (4.3) 

In App. A we prove that at least | Q\ of the p„'s are equal 
to unity. This immediately implies the inequality 

G/G >|e|. (4.4) 

For fe|| = we have, additionally, 

G/G > Co, for fc|| - 0. (4.5) 

In a topologically trivial system, with Q, Qo — 0, these 
inequalities are ineffective, while for \Q\, Qo — 2 the in- 
equalities are saturated (since G cannot become larger 
than 2Gq). Scattering events in the normal or super- 
conducting region that conserve feii, such as spin mixing, 
cannot change the conductance once it is saturated. 

V. EFFECTS OF ADDITIONAL UNITARY 
SYMMETRIES 

Further unitary symmetries may enforce restrictions 
on both the topological invariant and the angle-resolved 
conductance, or even introduce new topological invari- 
ants. In the first subsection we consider spatial sym- 
metries that invert feii —fey, whereas in the second 
subsection we address symmetries that conserve fe||. 

A. Spatial symmetries 

We consider a spatial symmetry of the form 

r(fe||) = (cr a <g> Tfe)r(-fe||)(<T a ® n). (5.1) 
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a, b 


X, X ox z,x 
or 0, y 


x, or 2,0 
or y, z 


x, z or y, x 
or z, z or y, 


y,y or 0, z 


x,y or z,y 
or 0, a; 


0,0 


Tab 


+ 1 


+ 1 


+ 1 


-1 


-1 


-1 


r 2 


+1 


+1 


-1 


-1 


+1 


+1 


class 


BDI 


BDI 


CI 


CII 


Dili 


Dili 


Qab 


0,±1,±2 


0,±1,±2 





0,±2 


0,2 


0,2 


G/Go 


> |Qo6| 


= |Qafc| 


X 


= |Qab| 


X 


= \Qab\ 



TABLE I: The first row lists the spatial symmetry (5.1); the second and third rows give the square of the anti- unitary operators 
(5.2) and (5.3); the fourth and fifth rows show the corresponding symmetry class and the values taken by the topological 
invariant; finally, the last row gives the relation between conductance and invariant for a topologically nontrivial system (so 
for Qab 7^ 0, with x indicating the absence of a relation). 



Combined with time-reversal symmetry (2.3) and charge- 
conjugation symmetry (2.4), this produces the two sym- 
metry relations 

r(fc||) =TabrHk ll )T ab \ T ab = [a a ■ (T y ) <g) Tfc/C, (5.2) 
r{k\\) = C ab r{k\\)C-^, C ab = a a ® (r b ■ t x )K, (5.3) 

where K is the operator of complex conjugation. The 
product of Tab and C a b brings us back to the chiral sym- 
metry (2.5). 

1. Topological invariant 

Depending on whether the anti-unitary operators T ab 
and C a b square to +1 or —1, the reflection matrix falls in 
one of the four Altland-Zirnbauer symmetry classes BDI, 
CI, CII, DHL 20 The various cases are listed in Table I. 
These all have a higher symmetry than the class AIII 
from which we started (with only chiral symmetry) . The 
additional symmetry may restrict the topological invari- 
ant to a smaller range of values. In class Dili a new Z 2 
topological invariant appears, that can be nonzero even 
if the Z invariant vanishes. 

We denote the modified topological invariant by 
Qab{k\\). In class CI only topologically trivial systems 
exist, 3 meaning that the spatial symmetry allows only 
for Qab = 0. For the other three symmetry classes the 
topological invariants are given by 14 

Q ab = ±Tri?e {-2,-1,0,1,2}, for BDI, (5.4) 
Qab = ±Tri?e {-2,0,2}, for CII, (5.5) 
Q ab = l + Pf(cr a ®T b ){a y r)e{0,2}, for DHL (5.6) 

The restriction to even integers in class CII (a 2Z in- 
variant) is a consequence of the Kramers degeneracy of 
the eigenvalues of the Hcrmitian matrix R = (a y ® r x )r. 
Symmetry class Dili has a Z 2 invariant. 



2. Conductance 

The expressions (5.4) and (5.5) for Q ab in class BDI 
and CII are the same as the expression (2.8) for Q in 



class AIII, so the topological invariant still provides a 
lower bound on the angle-resolved conductance, 



G/G >\Qab\, for BDI and CII. 



(5.7) 



In symmetry class Dili the invariants Qoo in Eq. (5.6) 
and Qo in Eq. (2.11) also have the same expression, so 
the inequality (4.5) still applies, 



G/G > Qoo- 



(5.8) 



No relation with the conductance exists for the other 
invariants in class Dili, so Qo x , Qxy, and Q zy provide no 
restriction on the conductance. 21 

The inequality (5.7) can be sharpened further in class 
BDI, so that it becomes an equality not only for | Q a b | = 2 
but also for |Q a &| = I- 22 As we show in App. C, this 
equality is enforced by the spatial symmetry (5.1) for 
(a, b) e {(y, z), (x, 0), (z, 0)}, so for three out of the six 
symmetries in class BDI. 

The last row of Table I summarizes the relation be- 
tween the topological invariant and the conductance for 
a topologically nontrivial system (Q a b 7^ 0). It is an 
equality for all symmetries in class CII and for some sym- 
metries in classes BDI and DHL 



B. Symmetries that preserve fey 

A different type of unitary symmetry preserves parallel 
momentum, 



r(fc||) = (a a n)r(k\\)(a a <g> n). 



(5.9) 



Combined with the chiral symmetry relation (2.5) and 
unitarity of r, this symmetry ensures that the matrix 
R = (<r a <£> Tb)R is a unitary matrix that squares to ±1. 
We can thus define a new Z invariant 



Q(feii) 



±Trfl(fc||) if R 2 = 1, 
|iTrfi(fe||) if R 2 = -1. 



(5.10) 



In general, Q and Q are distinct, and in particular Q 
can be an even function of fen. The coexistence of two 
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FIG. 2: Interface between a normal metal and a 2D Rashba 
superconductor. The Fermi surface is split into two circles, 
which intersect the nodal lines (red) of the superconducting 
pair potential in eight nodal points. 
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distinct topological invariants is quite unusual, and as 
we will see, it has observable consequences in the con- 
ductance. 

For b € {0, z} nonzero values of Q constrain the con- 
ductance in the same way that Q does in Eq. (4.4). For 
b G {x, y} one has instead the constraint 

G/G <2-\Q\, (5.11) 
as we show in App. B. 



VI. APPLICATION: 2D RASHBA 
SUPERCONDUCTOR 

As a first application of our general scattering theory 
we consider a two-dimensional superconductor with spin- 
singlet and spin-triplet pairing mixed by Rashba spin- 
orbit coupling. The topologically protected edge states 
for this Rashba superconductor have been studied in 
Refs. 10,23,24 using the Hamiltonian formulation. We 
summarize those results in the next subsection, before 
proceeding to the scattering formulation and the calcu- 
lation of the conductance. 



A. Hamiltonian and edge states 

The superconductor has the Bogoliubov-de Gennes 
Hamiltonian 



H(k) 



e(k)+g(k) 
k At(fe) 



A(k) 
-e(k)+g(k) 



with free electron part e(fc) = \k\ 2 /2m — /i, at Fermi 
energy /i, and Rashba spin-orbit coupling g(k) = 
X(k y ,—k x ,0). We have set ft = 1 and have collected 
the three Pauli matrices in a vector er = (a x ,a y ,a z ). 
The Fermi surface consists of two concentric circles at 
momenta 

k± = [{m\f + 2m(j] 1/2 ± mX. (6.2) 

For later use we give the spin orbit energy E so = mX 2 and 
the spin-orbit momentum and length fc so = raA = l/l so - 



FIG. 3: Topological invariant Q = -W of the 2D Rashba 
superconductor (<f> — 0, /i = 10 E so ), as a function of momen- 
tum k y along the NS interface and ratio At / A s of triplet and 
singlet pairing energies. 



The mixed singlet-triplet pair potential is given by 
g(k) ■ cr 



A(fc) = /(fc) A s + A 



X(2mii) 1 / 2 



1 



f(k) = 2^ [k x k y cos 2cj> + \(k 



k~ . 



(6.3) 
sin 20] , (6.4) 



The strength of the singlet and triplet pairing is param- 
eterized by the energies A s and A t . The nodal lines of 
vanishing pair potential are oriented at an angle <p with 
the NS interface (see Fig. 2). The intersection of the 
nodal lines with the Fermi surface defines 8 nodal points, 
in each of which Det H = 0. 
The chiral symmetry 

H(k) = -(<7 y ®T x )H(k)(<r y ®T x ) (6.5) 

ensures that H can be brought in the off-diagonal form 



U^H{k)U 



q(k) 
ot(fc) t 



(6.6) 



The Z topological invariant is then defined by the wind- 
ing number 10 

W(k y ) = i-Im J dk x -^lxLDetq(k x ,k y ), (6.7) 

for any k y that is not equal to the projection of one of 
the nodal points on the y-axis. 

As analysed in Refs. 10,23,24, the termination of the 
superconductor at x — by an insulator (or by vacuum) 
produces IW^fcy)! dispersionless edge states (flat bands). 
A simple example occurs for (f> — and A t = 0, corre- 
sponding to dry-wave spin-singlet pairing. Then 



W(ky) 



2sign(fc w ) if \k y \ < fc_, 



sign(fc y ) 




if fc_ < \ky\ < k 

if \k y \ > k + , 



+ • 



(6.8) 



so there are two topologically protected edge states for 
\k y \ < k- and a single one for k- < \k y \ < k + . 
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FIG. 4: Topological invariant Q of the reflection matrix from 
the 2D Rashba superconductor, as a function of momentum k y 
along the NS interface and angle 4> between the interface and 
the nodal line. The left panel shows results for spin-singlet 
pairing (A s = E so , At = 0) and the right panel for spin-triplet 
pairing (A t = E so , A s — 0). In both panels fi = 10i? a o and 
/in = 30 -Eso- The dotted lines indicate a topologically trivial 
system in class CI, as a consequence of the spatial symmetry 
(6.9). 



For nonzero A t the phase boundaries 
affected in the interval 



remain un- 



-y / 2mn/k- < A t /A s < y/2mfj,/k+, 

see Fig. 3. To contrast the spin-singlet and spin-triplet 
dominated regimes, we will in what follows focus on the 
two limits A t — > and A s — > 0. 




FIG. 5: Electrical conductance and Z topological invariant 
for three of the angles <j> from Fig. 4. A nonzero Z2 invariant 
appears in the spin-triplet case: Qo = 2 for k y = 0, 4> 7^ 0. 



B. Reflection matrix and conductance 

If the superconductor is not terminated at x — but 
connected to a normal metal, the edge states hybridize 
with the continuum of the metallic bands. The topolog- 
ical signature then shows up in the conductance rather 
than in the density of states. To reveal these signatures 
we construct the reflection matrix of the NS interface 
and calculate both the topological invariant (2.8) and 
the angle-resolved conductance (4.1). 

We used either an analytical method of calculation 
(matching wave functions at the NS interface), or a nu- 
merical method (discretizing the Hamiltonian (6.1) on a 
square lattice and calculating the Green function). We 
made sure that the lattice constant was sufficiently small 
that the two methods gave equivalent results. In the nor- 
mal metal we set both the pair potential and the spin- 
orbit coupling to zero, so that there is a single Fermi 
circle with momentum kjss — (2m /i-^) 1 / 2 . Because of a 
potential step at the NS interface, the chemical potential 
/in in the normal metal (x < 0) can differ from the value 
H in the superconductor (x > 0). 

Results are collected in Figs. 4 and 5. As a first check, 
we note that for (j> = 0, At = 0, we recover Eq. (6.8) - 
up to an irrelevant minus sign, Q = — W. For <fi — (n + 



l/2)7r/2, the system is topologically trivial, Q{k y ) = 0, 
regardless of the choice of A s , A t (black dotted lines in 
Figs. 4 and 5). This can be understood as a consequence 
of spatial symmetry: For cos 2<f> = the system fulfills 

H[k x ,k y ) = (jyH(k x , -k y )a y => r(k y ) = cr y r(~k y )<7 y . 

(6.9) 

This is a symmetry condition of the type (5.1), with 
a = y,b = 0, forcing the reflection matrix into the topo- 
logically trivial symmetry class CI (see Table I). At 
k y = the Z invariant Q vanishes, but the Z2 invari- 
ant Qo can be nonzero. This happens for A s = 0, <p 
(mod 7r/2), when Qo = 2. 

Fig. 5 shows how the topological invariant enforces the 
quantization of the angle-resolved conductance. First of 
all, G/G = 2 whenever \Q\ = 2 or Q = 2. For = 
quantized plateaus at G/Gq — 1 appear because of the 
spatial symmetry 

r(ky) = (<j y ® T z )r(~k y )((j y ® t 2 ), (6.10) 

which is a symmetry of the type (5.1) with a,b — y,z. 
This forces the reflection matrix into class BDI and en- 
sures that the conductance is quantized for any nonzero 
Q (see Table I). 
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FIG. 6: Topological invariants Q (left panel) and Q (right 
panel) for an NS junction between a normal metal and the 
anisotropic Rashba superconductor of Sec. VI C. The param- 



eters chosen are: A s = 
The Z2 invariant Qoo 
panel. 



5 so , At = 0, n = 10-Eso, a*n = 30 -E so . 
2 on the dotted red lines in the left 



C. Anisotropic spin-orbit coupling 

A strongly anisotropic dispersion, m x 3> m y , can 
produce an anisotropic spin-orbit coupling term of the 
form 25 g(k) = A(0,— k x ,0). Topological invariants and 
conductance are plotted for the spin-singlet regime (A t = 
0) in Figs. 6 and 7. There are two qualitative differences 
with the isotropic case of the previous subsections. 

First of all, for </> = mr/2 the regions with |<2(fcj,)| = 1 
are missing. This can be explained by the spatial sym- 
metry 



r(k y ) = T z r(-k y )r z , 



.11) 



of the type (5.1) with a, b = 0, z. As a consequence, see 
Table I, the topological invariant Q(k y ) becomes a 2Z 
invariant of class CII, excluding |<2(fcj,)| = 1. 

Secondly, there is a unitary symmetry a y r(ky)<7 y = 
r(k y ) that holds for all (p. This allows us to define an 
additional topological invariant, 



Q = ^Tra y R = iTr r x r, 



3.12) 



following Sec. VB. The topological invariants Q and Q 
are independent, in particular, Q(k y ) = Q(—k y ) while 
Q(ky) = —Q(—ky)- Each topological invariant Q and Q 
gives a lower bound on the conductance. This explains 
the diamond-shaped regions in the phase diagram with a 
quantized conductance G/Gq =2, enforced by |Q| =2. 

There is a third invariant: At <f> — (n + l/2)7r/2 the 
spatial symmetry r(k y ) = r(—k y ) places the reflection 
matrix in symmetry class DHL According to Eq. (5.6), 
the corresponding Z2 invariant Qoo = 2 on the dotted 
red lines in the phase diagram. 

This third invariant does not lead to additional con- 
straints on the conductance, since we already have Q = 2 
when Qoo = 2. But the two invariants Q and Q are both 
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FIG. 7: Electrical conductance and Z topological invariants 
for three of the angles <f> from Fig. 6. 



needed to explain the quantized conductance. The coex- 
istence of two topological invariants is an unusual feature 
of this system. 



VII. EFFECTS OF ANGULAR AVERAGING 
AND DISORDER 

It may be possible to measure the angle-resolved con- 
ductance G(ku ), 26 but one typically measures the angular 
average. Moreover, disorder is detrimental for the con- 
ductance quantization if it mixes parallel momenta with 
different values of the topological invariant. In this sec- 
tion we investigate whether signatures of the conductance 
quantization can survive the effects of angular averaging 
and disorder. 

We focus on the 2D Rashba superconductor of Sec. VI, 
for At = 0, (f) = 0, when the topological invariant is given 
by Eq. (6.8). The angular average of the conductance for 
an interface of width W is given by 



W f N 
Gns = TT I dk y G(ky). 

271 J-fcN 



(7.1) 



The reflection matrix, which determines G(k y ) via Eq. 
(4.1), is calculated numerically using the square lattice 
discretization of the Hamiltonian (6.1) (lattice constant 
a = 0.2 Z so , W = 32 Z so ) . Disorder is added to a strip 
— L < x < (L = 31.6 l so ) of the normal region by means 
of a random on-site potential, distributed uniformly in 
(— Uo/2, Uo/2). Results are averaged over 100 disorder 
realizations. 

In Fig. 8 we show the dependence of Gns on the Fermi 
momentum kj$ in the normal region. This is relevant 
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FIG. 8: Average conductance (7.1) of the NS junction as a 
function of the Fermi momentum &n in the normal region, for 
various disorder strengths. The 2D Rashba superconductor 
has a d^-wave pair potential (<j> — 0, A t = 0, A s = E so , 
H — 10 -Eso). Disorder strengths from top to bottom curve: 
U /E BO = 0,1,2,3,4,5. 
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FIG. 9: Differential conductance of the NS junction for vari- 
ous disorder strengths. The parameters for the superconduc- 
tor are the same as in Fig. 8. In the normal region we have 
fixed = 25 -Eso- Disorder strengths from top to bottom 
curve: U /E BO = 0, 2.5, 5, 7.5, 10. 



VIII. THREE-DIMENSIONAL 
SUPERCONDUCTORS 

A. Topological invariant for arc surface states 

The topological invariants considered so far, and the 
resulting constraints on the angle-resolved conductance, 
apply both to 2D and 3D nodal superconductors. In this 
section we discuss features that are specific for 3D super- 
conductors. The topological invariant Q(fey) of Sec. II B 
then counts dispersionless surface states, pinned to zero 
energy (the Fermi level) in a 2D region of parallel mo- 
mentum fc|| = {k\,k 2 ). The boundary of this flat band 
region is formed by nodal rings, closed contours of fen on 
which transmission through the superconductor is possi- 
ble — in other words, the superconducting gap vanishes 
for k = fey). 

The new feature that appears in a 3D superconductor 
is the possibility of zero-energy boundary states along a 
ID arc connecting two nodal rings. Some aspects of their 
topological nature have been discussed in the Hamilto- 
nian formulation of Rcf. 19. Here we consider the alter- 
native scattering formulation, and use it to obtain topo- 
logical constraints on the conductance. 

We consider a spatial symmetry on the 2D surface of 
a 3D superconductor, in which only one of the two com- 
ponents of parallel momentum is inverted: 

r(k 1 ,k 2 ) = (<J a (g) T b )r(-k 1 ,k 2 )(cra ® n). (8.1) 

Along the line k 2 — 0, this is a symmetry of the type 
(5.1), so we can follow Sec. VAl to introduce topologi- 
cal invariants Q a b{k\). The resulting constraints on the 
angle-resolved conductance G(/ci,0) are summarized in 
Table I. 

Alternatively, for k\ — 0, the symmetry (8.1) is of 
the type (5.9) with topological invariant Q(k 2 ) from Eq. 
(5.10). The corresponding constraints on the conduc- 
tance are discussed in Sec. VB. 



if the normal region is a semiconductor, where one can 
vary &n by a gate voltage. The quantization of G(k y ) 
manifests itself as a quantized slope of Gns versus k^: 
the steep slope for ku < k- (where \Q\ =2) is reduced 
by a factor of two in the interval fc_ < &n < k + (where 
\Q\ = 1), and then is strongly suppressed for &n > k + . 
This signature of the topological invariant gradually dis- 
appears with increasing disorder. 

Another signature can be seen for fixed /cn in the de- 
pendence of the differential conductance dl/dV on the 
applied voltage V. As shown in Fig. 9, the peak in dl/dV 
around V — is a superposition of two peaks with dif- 
ferent widths, the narrower one originating from parallel 
momenta in the |Q| = 2 regions and the broader one 
from the \Q\ = 1 regions. The single edge state of the 
latter regions couples more strongly to the continuum of 
the metal and thus has a larger width. 



B. Example 

As an example, we apply these general considerations 
to the same Rashba Hamiltonian (6.1), but now with a 
3D dispersion, 

e{k) = {kl + kl + k 2 z )/2m- l i. (8.2) 

In the pair potential (6.3) we set f(k) = 1. This Hamil- 
tonian applies to non-centrosymmetric s+p-wave super- 
conductors of point group Ca v . As described in Ref. 19, 
these superconductors have arc surface states connecting 
two nodal rings. They appear for example for the (011) 
surface orientation that we will consider in the follow- 
ing. The two components of parallel momentum on the 
surface are k\ = k x and k 2 — (k y — k z )/y2. 

We can obtain two topological invariants from the re- 
flection matrix r(ki, k 2 ), plotted in the left panel of Fig. 
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FIG. 10: Topological invariants and electrical conductance 
of an NS junction between a normal metal and an s+p-wave 
superconductor of point group d v . The invariants Q and 
Q are plotted in the left panel. The right panel shows the 
conductance (black curve) and invariant Q (red dotted curve) 
along the line fci = 0. Parameters chosen are A s = 0.1 E so , 
A t = 0.2 E BO , n = 10E SO and ju N = 30£ so . 

10. The first invariant 

Q{kxM) = fca) = 5Tr( ( j y ®T 2: )r(fc 1 ,fc 2 ) (8.3) 

follows from chiral symmetry, see Sec. II, and is defined 
on the entire 2D plane of parallel momenta. This Z 
invariant is nonzero inside the regions bounded by the 
nodal rings, where it identifies a surface flat band. 

A second Z invariant appears as a consequence of the 
spatial symmetry 

H(k x ,k y ,k z ) = (<r x ® r z )H(-k x , k v , k z )(a x ® t z ) => 
r(fci,fc 2 ) = {a x ®T z )r(-k 1 ,k 2 )((T x ®T z ). (8.4) 

The line k\ — connects the two nodal rings and on this 
line the invariant 

Q{k 2 ) = ±Tr(<r x ®T z )R(0,k 2 ) = -iTr(a z ®r a )r(0,fc 2 ) 

(8.5) 

can take on a nonzero value. 

The non-trivial invariants enforce a lower bound on the 
conductance, as is illustrated in the right panel of Fig. 
10. This leads to a quantized conductance G/Gq = 2 
along the line k\ = 0. 

The symmetry (8.4) produces arc surface states on 
all surfaces parallel to the ^-direction. For the (010) 
surface analyzed in Ref. 19 there is an additional spa- 
tial symmetry, r{k Xl k z ) — r(k x ,—k z ). For k x = 
this additional symmetry allows for the Z 2 invariant 
Qoo = 1 + Pf <7yr(0, k z ), in addition to the Z invariant 
(8.5). For other surface orientations (Onm) only the Z 
invariant is responsible for the arc states. 

IX. CONCLUSION 

In conclusion, we have constructed a topological in- 
variant Q(fc||) of the Andreev reflection matrix at the 
interface between a time-reversal symmetric nodal su- 
perconductor and a normal metal. In the absence of a 



tunnel barrier, this interface has no zero-energy bound- 
ary states, but the topologically nontrivial phase can still 
be detected in the angle-resolved conductance G(ku). A 
variety of symmetry classes can be realized (AIII, BDI, 
CI, CII, Dili), by allowing for additional unitary symme- 
tries. The corresponding topological invariants are given 
by a trace or Pfaffian of the reflection matrix. 

Many of these topological invariants have been stud- 
ied before in the Hamiltonian formulation for an infinite 
system. 10,11,16-19 The scattering formulation presented 
here makes it possible to directly relate Q(feii) to G(kn). 
We have systematically examined when a nontrivial topo- 
logical invariant enforces a quantized conductance, and 
when it only provides a lower bound. This approach can 
identify surface flat bands (within nodal rings) as well as 
arc states (connecting nodal rings), even when these zero- 
energy boundary states have merged with the continuum 
of states in the normal metal. 

We have applied the general theory to 2D and 3D su- 
perconductors with spin-singlet and spin-triplet pairing 
mixed by Rashba spin-orbit coupling. The appearance of 
a quantized conductance has allowed us to verify known 
topological invariants and to identify new ones. In par- 
ticular, in the 2D case of a strongly anisotropic spin-orbit 
coupling, we have shown the coexistence of two topolog- 
ical invariants — which provide independent constraints 
on the conductance. 

To make contact with experiments, the effects of angu- 
lar averaging and impurity scattering on the conductance 
quantization have been investigated by numerical simu- 
lation of a disordered NS interface. 
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Appendix A: Topological invariant counts number of 
unit Andreev reflection eigenvalues 

1. Proof for the Z invariant 

The Hermitian matrix R% h has eigenvalues p n £ [0, 1]. 
We wish to prove that at least \Q\ of these Andreev re- 
flection eigenvalues are equal to unity. 

Let 4> be an eigenvector of R e t with eigenvalue A. As- 
sume A ^ ±1 (so p = A 2 < 1). Since R 2 eh = 1 - R} hh R hh , 
the vector <f>' = Rhh4> cannot vanish. Since RheRhh = 
—RhhRehi it then follows that <j)' is an eigenvector of Rhe 
with eigenvalue p = —A. 
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Now consider the Z topological invariant Q = 
\ Ent'' 1 " + A 4 ") m symmetry class AIII. The eigenvalues 
A„ 7^ ±1 of R e h are cancelled by an eigenvalue fi n = — A„ 
of i?/, e . The cancellation can only be avoided for the M 
eigenvalues A„ equal to ±1, resulting in \Q\ < M — as 
we set out to prove. 

2. Proof for the Z2 invariant 

For any 4x4 antisymmetric matrix A with a block 
structure, 




the Pfafhan is given by 

Pf A = -DetAis - ITyA u A 22 . (A2) 

We apply this identity to the antisymmetric matrix 
a y r at kii — 0, to obtain the Z 2 topological invariant in 
symmetry class Dili, 

Qo = 1 - Det R eh - \ Tr R ee R hh 

= l-DetR eh -±Tv(l-R 2 eh ). (A3) 

In terms of the two eigenvalues Ai,A2 G [ — 1, 1] of R e h 
this reduces to 

Qo-|(Ai-A 2 ) 2 . (A4) 

Since by construction Q equals either or 2, we have 
either Q Q = Ai = A 2 or Q Q = 2 ^ X 1 = -A 2 = ±1. 
This shows that at least Qo of the Andreev reflection 
eigenvalues p n = A 2 are equal to unity. 

Appendix B: Proof of Eq. (5.11) 

We consider the topological invariant (5.10), con- 
structed from the matrix R = (a a (&Tt,)R with b e {x, y}, 
and wish to proof the constraint (5.11) on the conduc- 
tance. This amounts to a proof that at least \Q\ of the 
Andreev reflection eigenvalues are equal to zero. 

We define the Hcrmitian matrix 

R= (Ree Reh] = { R(k\\) if i? 2 = 1, (m) 

\R he Rhh) \iR(k\\) if R 2 = -l. 



Let 4> be an eigenvector of R ee with eigenvalue A. Assume 
A 7^ ±1. Since R 2 e = 1 — W he Rh e , the vector <fr' = Rhefy 
cannot vanish. With Rh e R ee — —RhhRhe, it then follows 
that 4>' is an eigenvector of Rhh with eigenvalue fi = —A. 

Now since Q = \ r Tv{R ee + R hh ) = \ X) n (A„ + fJ. n ), the 
eigenvalues A„ 7^ ±1 of R ee are cancelled by eigenvalues 
jj, n = — A„ of Rhh in the expression for the topological 
invariant. The cancellation can only be avoided for the 
M eigenvalues A„ equal to ±1, resulting in | Q\ < M . The 
existence of at least \Q\ unit eigenvalues of R\ e R ee = R 2 ee 
is equivalent to the existence of at least | Q| zero Andreev 
reflection eigenvalues and thereby proves Eq. (5.11). 
Appendix C: Equality of conductance and 
topological invariant in class BDI 



A topologically nontrivial 4x4 reflection matrix in 
class BDI has either |Q a b| = 2 or |Q a &| = 1- In the 
former case the inequality (5.7) is saturated, because 
G/Gq < 2, but in the latter case it provides only a lower 
bound on the conductance. We now wish to show that 
the inequality can be sharpened to an equality for three 
of the six spatial symmetries (5.1) in class BDI. More 
precisely, we will show that \Q a b\ = 1 implies G/Go = 1 
for (a,b) e {(y,z),(x,0),(z,0)}. 

For each of these three cases the symmetry relation 
(5.2) implies that R he = cr a Rj h (j a , so Tr R he = Tt R eh - 
Denote the eigenvalues of R e h and Rhe by Ai,A 2 and 
/ii, /i 2 , respectively. (All are real numbers in the interval 
[—1, 1].) The equality of the traces gives Ai + A 2 = + 
/i 2 . The topological invariant (5.4) determines the sum 
Ai + A 2 + ^ + fi 2 = ZQab, hence Ai + A 2 = Q ab . 

Because classes BDI and AIII have the same expression 
for the topological invariant, we may apply the result of 
App. AI that at least |Q a fc| of the A„'s equal ±1. If we 
take \Q a b\ = I? I Ai| = 1, then necessarily A 2 = 0. The 
dimensionless conductance G/G = \\ + A 2 thus equals 
unity, as we set out to prove. 

Our finding can be seen in a broader context as a 
manifestation of Beri degeneracy of Andreev reflection 
eigenvalues: 27 The charge-conjugation symmetry (5.3), 
with (a, b) e {(y, z), (x, 0), (z, 0)}, enforces a twofold de- 
generacy of the Andreev reflection eigenvalues p n = A 2 
that can only be avoided if p n equals or 1 . 
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